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On the existence of a first order phase transition at small vacuum angle 9 
in the CP 3 model* 
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We examine the phase structure of the CP 3 model as a function of 9 in the weak coupling regime. It is shown 
that the model has a first order phase transition at small 9. We pay special attention to the extrapolation of the 
data to the infinite volume. It is found that the critical value of 9 decreases towards zero as fi is taken to infinity. 



1. INTRODUCTION 



die away like 1 / N [5] . 



One of the characteristic features oiQCD is the 
existence of an integer valued topological charge 
Q. This gives rise to an additional term in the 
action, 



Sg = S + i6 



G [0,2tt), 



(1) 



where each value of 6 results in a different vacuum 
state. A priori 6 is a free parameter. Since the 
vacuum is an eigenstate of CP only for 6 = 0, 7r 
and no CP violation is observed in the strong in- 
teractions, 6 must, however, be very close to zero. 

Is there any dynamical reason for that? There 
are indications that confinement at 6 = is 
caused by the condensation of color magnetic 
monopoles [1]. In the 6 vacuum these monopoles 
acquire an electric charge 6/2tt [2]. One may then 
argue [3] that a rich phase structure will emerge 
as a function of 6. In the weak coupling limit it 
is even conceivable that confinement will be lim- 
ited to 6 = only which would solve the problem 
rather elegantly. 

Since many of the characteristic features of 
QCD are shared by the CP N_1 models, it sug- 
gests itself to examine these models first. We have 
chosen the CP 3 model. For smaller N one may 
run into the problem of dislocations [4], while for 
large N the signal one hopes to see will probably 
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2. LATTICE SIMULATION 

We use the action 
S = -2f3j2 TrP ( x ) p ( x + A)> 



(2) 



where P(x) = z(x)z^(x), TrP(x) = |z(a;)| 2 = 1, 
and we employ periodic boundary conditions. 
The action can be rewritten in terms of the com- 
plex scalar field z(x) and a composite abelian 
gauge field interacting minimally with each other. 
The associated U(l) link field has the form 



(x) ■ z(x + fl) 
\z^(x) ■ z(x + fl)\ 



11 = 0,1. 



(3) 



The topological charge is obtained from (3) by 



forming plaquette fields U p 
7T and writing 



Q 



-7T < F 01 < 



(4) 



As can easily be seen, Q is an integer. 

In this talk we will restrict ourselves to the cal- 
culation of the partition function and a few quan- 
tities which derive from it. A full account of our 
work will be given elsewhere [6]. The partition 
function reads 



Z(ff) 



Q 



P(Q), 



(5) 



1 



p(Q) 



JJVzVz^QHjzl - l)e- 
JVzVzH(\z\ 2 - l)e-' 



(6) 



where the subscript Q indicates that the path 
integration is restricted to the given charge sec- 
tor. In order to compute p(Q) we proceed as fol- 
lows [7]. We divide the phase space into over- 
lapping sets of five consecutive charges. In each 
of these sets the z fields are updated by a com- 
bination of Metropolis and overrelaxation steps, 
z' = exTp(i(f) a X a )z, where the ?7 (4) generators A a 
are selected randomly [8]. This is supplemented 
by an additional acceptance criterion: if the new 
charge is in the same set the configuration is ac- 
cepted and the new charge recorded; if the new 
charge is outside the set the change is rejected 
and the old charge recorded. Furthermore, a trial 
charge distribution which is approximately equal 
to the true distribution is incorporated. 

So far we have done simulations at two values 
of [3 on various lattice volumes. We have chosen 
/3 = 2.5, V = 28 2 ,32 2 ,38 2 ,48 2 , and /? = 2.7, 
V = 46 2 ,56 2 ,64 2 ,72 2 . The correlation lengths 
at these two [3 values are [9] ^ ps 4.5 and ^ ps 
8.8, respectively. On these lattices we are able 
to follow p(Q) over more than twenty orders of 
magnitude. 

3. RESULTS 

The first quantity we looked at is the free en- 
ergy F(9) which is given by 



F(0) 




8/u 
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Figure 2. The average topological charge q{9). In 
the second half of the 9 interval not shown here 
q(6) = -q(2w-6). 



Z(0) 



-VF(6) 



(7) 



A first order phase transition will manifest itself 
in a kink in F(9). In Fig. 1 we show F(9) on the 
V = 64 2 lattice at (3 = 2.7. The error bars have 
been computed by a jackknife method. We see 
two kinks, one at 9 ps 0.57T (and a corresponding 
one at 9 ps 1.57r) and one at 9 = 7r. This should be 
compared with the prediction of the dilute instan- 
ton gas approximation, F(9) oc 1 — cos^, and the 
prediction of the l/N expansion [5], F(9) oc 9 2 . 
Clearly, our results do not agree with either of 
the two. In the following we shall mainly be in- 
terested in the phase transition at small 9. 

The first derivative of F(9) is equal to the av- 
erage topological charge q(9): 



dF(9) 
A9 



1(0) 



1 



V Z{9) 



J2^ sq Qp(Q)- (8) 

Q 



Figure 1. The free energy F(9). 



At a first order phase transition this quantity 
must be discontinuous. In Fig. 2 we show q(9) 
for the V = 64 2 lattice at j3 = 2.7. We see 
that q(9) increases almost linearly with 9 until 
it reaches the location of the phase transition 
where it jumps to a value which is close to zero. 
According to (4), q(9) can be interpreted as a 
background electric field. In Ref. [6] we shall 
argue that at the phase transition the charged 



2 
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Figure 3. The free energy F(9). The error bars 
are omitted in order to keep the figure legible. 
They are comparable to the ones in Fig. 1. 



Figure 4. The critical value of i 



4. CONCLUSIONS 



scalars (associated with the field z(x)) are liber- 
ated which causes the background electric field to 
collapse. 

Let us now look at the volume dependence of 
the free energy. In Fig. 3 we show F(9) for four 
different volumes at [3 = 2.7. We see that the 
location of the phase transition moves to smaller 
values of 9 as V is increased. Below the phase 
transition all data points fall on a universal curve. 
If the lattice volume is too small we see no signal 
of a phase transition [10]. 

We denote the critical value of 9 by 9 C and iden- 
tify it with the local minimum of Z(9). The vol- 
ume dependence of 9 C is linked to the order of the 
phase transition. For a first order transition we 
expect 

9 c (V)-9 c (x>)ocV- 1 . (9) 

In Fig. 4 we show 9 C as a function of V~ x for (3 
= 2.5 and 2.7. In both cases our data fall on a 
straight line in accordance with a first order phase 
transition. This allows us to extrapolate the data 
to the infinite volume. 

We obtain c (oo) = 0.32(2) at /? = 2.5 and 
6* c (oo) = 0.18(3) at [3 = 2.7. It is perhaps not 
surprising that the critical value of 9 depends 
strongly on [3. Note that 9 C = 7r in the strong 
coupling limit [11]. 



We have found a novel, first order phase tran- 
sition in 9 with 6' c (cxd) decreasing towards zero as 
we approach the continuum limit. This transition 
is consistent with a deconfining phase transition 
[6]- 

If somebody wants to repeat the calculation he 
should start with j3 = 2.7, V = 56 2 . 
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